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ABSTRACT 

In this paper we study the supercluster - cluster morphological properties 
using one of the largest (2 x 512 3 SPH+N-body simulations of large scale 
structure formation in a ACDM model, based on the publicly available code 
GADGET. We find that filamentary (prolate- like) shapes are the dominant 
supercluster and cluster dark matter halo morphological feature, in agreement 
with previous studies. However, the baryonic gas component of the clusters is 
predominantly spherical. We investigate the alignment between cluster halos 
(using either their DM or baryonic components) and their parent supercluster 
major-axis orientation, finding that clusters show such a preferential align- 
ment. Combining the shape and the alignment statistics, we also find that the 
amplitude of supercluster - cluster alignment increases although weakly with 
supercluster filamentarincss. 

Key words: galaxies: clusters: general - cosmology: theory - large-scale struc- 
ture of universe 



1 INTRODUCTION 

The wealth of recent observational data has brought 
great progress in understanding the cosmic structure 
formation pattern. At the largest cosmic scales it is 
known that clusters of galaxies are not randomly dis- 
tributed but tend to aggregate in even larger conglom- 
erations, the so called superclusters (cf. Oort 1983; Bah- 
call 1988). The large sizes of superclusters (J> 30 h~ l 
Mpc) in conjunction with the amplitude of their mem- 
ber cluster peculiar velocities (~ 1000 km/sec) imply 
that superclusters constitute unvirialized density fluc- 
tuations which should reflect the initial conditions that 
gave rise to structure formation processes. 

Observational and theoretical studies of the super- 
cluster morphology and environment suggest that they 
are not spherical but instead their shapes are elongated 
(eg. Zeldovich, Einasto & Shandarin 1982; de Lappar- 
ent, Geller & Huchra 1991) with a prolate- like ten- 
dency (eg. Frenk et al. 1988; Plionis, Valdarnini & Jing 
1992; Dubinski 1994; Einasto et al. 2001; Jing & Suto 
2002; Diaferio et al. 2003; Sheth et al. 2003; Einasto 
et al. 2003; Shandarin, Sheth, Sahni 2004; Pandey & 



Bharadwaj 2005). Sathyaprakash, Sahni & Shandarin 
(1998), Basilakos, Plionis & Rowan-Robinson (2001), 
Kolokotronis, Basilakos & Plionis (2002), using either 
the IRAS galaxy distribution or the Abell/ACO clus- 
ters, gave attention to the cosmological inferences of su- 
percluster shape statistics claiming that a low matter 
density flat cosmological model (f2 m = 1 — Qa = 0.3) 
fits the large-scale observational results at a high signif- 
icance level. 

A variety of indications support the formation of 
clusters by hierarchical aggregation of smaller units 
along filamentary large-scale structures (eg. West 1994; 
Ostriker & Cen 1996 and references therein). Under this 
scenario dynamically young clusters will have a ten- 
dency to be aligned with neighboring structures, since 
the accretion of matter takes place along the large- 
scale filamentary structure within which the clusters 
form. Indeed, observational studies of structure orien- 
tations, which has a long history in cosmology, have 
found strong indications of various alignment effects. 
Bingelli (1982) was the first to find that the major axes 
of neighboring clusters of galaxies tend to point toward 
each other. On the other hand, Struble & Peebles (1985) 
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have failed to measure a significant alignment signal (see 
also Flin f987; Rhee & Katgert 1987; Ulmer, McMillan, 
& Kowalski 1989). However, in the last decade many 
authors have claimed that cluster formation processes 
are strongly connected to the supercluster network and 
thus generate measurable environmental effects, among 
which strong alignment effects observed in both obser- 
vational and N-body data (eg. West 1994; Plionis 1994; 
West, Jones, & Forman 1995; van Haarlem, Frenk & 
White 1997; Colberg et al. 2000; Onuora & Thomas 
2000; Chambers, Mellot & Miller 2002; Plionis & Basi- 
lakos 2002; Faltenbacher et al. 2002; Plionis et al. 2003; 
Bailin & Steinmetz 2004; Kasun & Evrard 2004; Pimb- 
blet 2005; Hopkins, Bahcall & Bode 2005; Lee, Kang & 
Jing 2005; Faltenbacher et al. 2005). 

Furthermore, Plionis (2002; 2004) using the APM 
galaxy data, has shown that the alignment effects are 
not confined only between cluster pairs or between lu- 
minous galaxies and their parent clusters, but extend to 
alignments between clusters and their parent superclus- 
ter major axis. This again is the sort of picture that one 
may expect in hierarchical structure formation scenar- 
ios where matter and galaxies flow within one or two 
dimensional structures, on the intersects of which the 
clusters form. 

The aim of this work is along the same lines, at- 
tempting to make a detailed investigation of the con- 
nection between the morphological and environmental 
properties of the large scale structure in the concor- 
dance ACDM cosmology. For a detailed study of the 
morphological properties of the supercluster-void net- 
work in such a cosmological model we refer the reader 
to the recent work of Shandarin, Sheth & Sahni (2004). 

The plan of the paper is as follows. The simulated 
cluster samples are presented in section 2. In section 
3, we briefly describe the method used to investigate 
supercluster and cluster shape properties and comment 
on some systematic effects related to the definition of 
superclusters. In section 4, we discuss the large scale 
structure orientation effects and finally in section 5, we 
present our conclusions. 



2 THE SIMULATED ACDM CLUSTERS 

In this study we use large-scale N-body simulations 
of the popular ACDM cosmological model in order to 
quantify the 3D cluster and supercluster morphological 
properties and their relation to the large scale structure 
network. 

2.1 N-body Simulation 

We have performed numerical simulations of a 500 ft -1 
Mpc cubic volume in which a random realisation of 
the concordance ACDM (fi m = 0.3, fl A = 0.7, h = 
0.7, Q,b = 0.045, as = 0.9) power spectra was generated 
with 2048 3 particles. Due to computational limitations, 
we have not yet simulated down to these resolutions but 
only selected regions. For the whole box, we increased 
the mass resolution in steps of a factor of 8 by averaging 
2048 3 initial distribution of particles. Then we replace 



each particle by a dark matter and gas particle. We have 
run simulations with 2 x 256 3 and 512 3 particles. Since 
all simulations were extracted from the same initial con- 
ditions, the same structures are formed, so that we can 
estimate the effects of mass resolution on the results (see 
§3.2.1). 

The simulations were done with an updated ver- 
sion of the parallel Tree-SPH code GADGET (Springel, 
Yoshida & White 2000). The code uses an entropy- 
conserving formulation of SPH (Springel & Hernquist 
2002) which alleviates problems due to numerical over- 
cooling. In order to study the gas dynamics of clusters, 
we ran the simulations with the same number of sph 
and dark matter particles. The results reported in this 
paper are based on the highest resolution simulation car- 
ried out which consists of 2 x 512 3 (~ 10 s ' 43 ) particles 
with a mass resolution of rridark = 6.6 x 10 10 ft _1 M Q and 
m sph = 1.2 x 10 10 /i _1 M Q . We still have a factor of 64 in 
mass left until we reach the resolution limit of the initial 
conditions. This allows further improvement in numeri- 
cal results on larger computers. Anyway, our simulation 
is already one of the largest adiabatic SPH simulations 
of large scale structure done so far. The mass resolution 
of the 2 x 512 3 simulation allows us to reliably iden- 
tify from large galactic halos (100+100 particles) to the 
biggest galaxy clusters (4 x 10 4 particles). 

The spatial force resolution was set to an equivalent 
Plummer gravitational softening of 15 h^ 1 comoving 
kpc. The sph smoothing length was set to the distance 
to the 40 th nearest neighbour of each sph particle. In 
any case, we do not allow smoothing scales to be smaller 
than the gravitational softening of the gas particles. 

These experiments have also been used in order to 
re-simulate individual clusters found in the low resolu- 
tion run with full resolution to study the properties of 
the ICM in hot X-ray clusters (Yepes et al. 2004). The 
number density of clusters found in this simulation agree 
quite well with recent observations of the X-ray temper- 
ature function by Ikebe et al. (2002) ranging from 1 keV 
to 10 keV clusters (for further details see Yepes et al. 
2004). 

The simulations were performed in the IBM Re- 
gatta p690+ supercomputer cluster at the John von 
Neumann Center Jiilich. For the largest simulation of 
2 x 512 3 particles, we employed 64 cpus simultaneously 
for a total of 64 x 17 = 10000 CPU hours. 

2.2 Simulated Cluster Samples 

At first we have calculated the minimal spanning tree 
(MST) of the DM particle distribution. The minimal 
spanning tree of any point distribution is a unique, well 
defined quantity which describes the clustering prop- 
erties of the point process completely (e.g., Bhavsar 
& Splinter 1996, and references therein). The minimal 
spanning tree of N points contains N — 1 connections. 
We are using an MPI program which calculates the MST 
of the 512 particles using 8 CPUs within about 10 min- 
utes. In a second step we cut the MST using different 
linking lengths in order to extract catalogs of friends- 
of-friends particle clusters. Note, that cutting a given 
MST is a very fast algorithm. We start with a linking 
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Figure 1. The mass function of FOF objects in the 500fe _1 
Mpc box with 512 3 dark matter particles (solid line) and the 
Shcth-Tormen approximation (dotted line) 

length of 0.17 times the mean inter particle distance 
which corresponds roughly to objects with the virializa- 
tion overdensity p/pmcan — 330 at [z = 0). Decreasing 
the linking length by a factor of 2 n (n =1,2,...) we get 
samples of objects with roughly 8" times larger over- 
densities which correspond to the inner part of the ob- 
jects of the first sample. With this hierarchical friends- 
of-friends algorithm we can also detect substructures of 
clusters (cf. Klypin et al. 1999). 

With a linking length of 0.17 we find ~ 77000 
objects in the simulation box with more than 100 
dark matter particles , i.e. total masses m fo f > 8 x 
10 12 ft _1 Mq. The corresponding mass function is shown 
in Fig. 1 together with the Sheth and Tormen (1999) 
analytical approximation. 

We proceed in the same manner to study the distri- 
bution of gas particles. In particular, we run the MST 
procedure over all gas particles to obtain catalogs of gas 
objects. Then we identify how gas and DM clusters are 
related. In most cases this is straightforward, i.e. the 
centers coincide. In a few cases (~ 2%) centers do not 
coincide. The simple reason is that the friends-of-friends 
algorithm may connect the more smoothly distributed 
gas particles by bridges which is not necessarily the case 
for the corresponding DM particles. 

For our present analysis we define two cluster 
samples based on two mass thresholds: (a) M c \ > 
1.34 x 1O 14 /i~ 1 M (hereafter ACDMi sample) and (b) 
M c i > 6.2 x lO 13 ft _1 M (hereafter ACDM 2 sample). 
These two subsamples contains 2773 and 7869 clus- 
ter entries with corresponding mean densities of ni ~ 
2.1 x 10~ 5 ft 3 Mpc~ 3 and n 2 ~ 6.1 x 10~ 5 ft 3 Mpc" 3 , re- 
spectively. 

These mass thresholds were chosen in order to 
roughly reproduce the spatial density of the REFLEX 
(Bohringer et al. 2001) and APM cluster samples (Dai- 




Figure 2. The spatial two-point correlation function of the 
ACDMi (M lim = 1.34 x lO 14 r 1 M ) and ACDM 2 (M Um = 
6.2 x 10 1[i h~ 1 Mq) samples. The error bars are estimated 
using the bootstrap procedure. The dashed lines represent 
the best-fitting power low £(r) = (r /r) 7 (see parameters in 
Table 1). 



ton et al. 1994). Utilizing the classical correlation esti- 
mator described by Efstathiou et al. (1991), we evalu- 
ate the real-space two-point correlation function, £(r), 
in logarithmic intervals. The resulting correlation func- 
tion parameters, for both cluster samples, are presented 
in Table 1. The derived correlation function slope is 
very near to its nominal value of 7 = 1.8 and we find 
r = 15.3 ± 1.4/i _1 Mpc and r = 11.8 ± 1.0/i _1 Mpc, 
respectively. It is clear that the correlation length in- 
creases with cluster richness, as expected from the well- 
known richness dependence of the correlation strength 
(eg. Bahcall & Burgett 1986; Bahcall & West 1992). In 
Fig. 2 we plot the derived two-point correlation function 
as filled dots, while the best-fitting power law model 
£(r) = (r /r) 7 is shown as straight lines (see Table 1). 

The ACDMi sample has a correlation length that 
within ~ 1.5a resembles that of the REFLEX X-ray 
cluster sample (Collins et al. 2000) while the ACDM 2 
sample has a correlation length which approaches that 
of the APM clusters (Dalton et al. 1994) and of a poor 
subsample of the SDSS-CE clusters (Goto et al. 2002) 
analysed by Basilakos & Plionis (2004). 



3 STRUCTURE SHAPE DETERMINATION 

3.1 Defining ACDM Superclusters 

The identification of superclusters is based on the use 
of the friend of friends algorithm (otherwise called per- 
colation algorithm; see Zeldovich et al. 1982), applied 
on the distribution of simulated clusters. The algorithm 
starts by placing a sphere, of a certain radius, around 
each cluster and then connects all neighboring spheres 
having an overlap region. Doing so for all clusters in the 
sample the algorithm provides a unique, for each specific 
percolation radius, list of connected clusters, dubbed 
"superclusters". It is evident that different percolation 
radii will result in different catalogues of superclusters. 
Therefore the choice of an optimal percolation radius 
is essential for the detection of superclusters which are 



© 0000 RAS, MNRAS 000, 000-000 



4 Basilakos et al. 



Table 1. List of the major parameters in the two analysed samples. The first column gives the model, column 2 gives the mass 
threshold, columns 3 and 4 give the clustering properties, column 5 gives the number of the simulated clusters contained in the 
volume used. Column 6 gives the number of the superclusters with 10 or more members (N^ p ) and finally, column 7 gives the 
percolation radius used. Note that r a , R pI have units of /i -1 Mpc and M lim has units of h~ 1 MQ. 



sample 


M lim /h- 1 M Q 


r 


7 


N c i 


N 10 

sup 


Rpr 


ACDMi 


1.34X10 14 


15.3±1.5 


1.82±0.07 


2773 


22 


18.0 


ACDM 2 


6.2xl0 13 


11.8±1.0 


1.78±0.04 


7869 


54 


11.5 




Figure 3. The 3D distribution of our detected superclusters 
(based on the ACDMi cluster sample). A Gaussian smooth- 
ing has been applied with lcr radius of 6 h^ 1 Mpc. The 
elongated prolate-like nature of the detected superclusters 
is evident. 



related uniquely to the specific underline point (cluster) 
process. 

To this end, we use a "critical" value of the per- 
colation radius (R pT ), which is related uniquely to the 
clustering properties of the underline cluster distribu- 
tion, given by Peebles (2001): 



Rt, 



'7 



(1) 



where uj s is the solid angle covered by the cluster sample, 
n is the mean cluster number density and r , 7 are the 
clustering parameters of the cluster sample (see Table 

!)■ 

Thus, for the ACDMi and ACDM2 cluster sam- 
ples, the critical percolation radius is R pT ~ 18 and 
11.5 h~ Mpc, respectively. Using these values, we ap- 
ply the percolation algorithm to the whole simulation 
box distribution of clusters, while in order to minimize 
edge effects, related to the simulation box boundaries, 
we rap-around whenever a supercluster touches the sim- 
ulation borders. The resulting number of superclusters, 
with 10 or more cluster members, is 22 and 54, respec- 
tively for the ACDMi and ACDM2 cluster samples. In 
Fig. 3 we present a three dimensional view of the de- 
tected (ACDMi) superclusters after smoothing them 
with a Gaussian window of a 6 h^ 1 Mpc radius. 

We have performed a further test in order to verify 
whether the detected superclusters are the optimal ones. 




15 20 
R pr (h-'Mpc) 

Figure 4. The probability that the supercluster detected in 
the simulated data with 10 or more members is real as a 
function of the percolation radius R pr . The solid line is for 
the ACDMi sample, the dashed for ACDM2. 



To this end we compare the number of detected ACDM 
superclusters for different values of the percolation ra- 
dius with those detected, for the same radius, in random 
samples of clusters having the same space number den- 
sity as the simulation clusters. We run a large number of 
Monte-Carlo simulations in which we destroy the intrin- 
sic ACDM clustering by randomizing the spatial coordi- 
nates of the clusters. On this intrinsically random clus- 
ter distribution, we apply the procedure described be- 
fore and identify the random superclusters, AT rail , which 
are due to our supercluster-identification method itself. 

We define the probability that the detected su- 
percluster in the simulated data is real as V = 1 — 
A?ran/A? SU p. In Fig. 4 we plot that probability depending 
on the percolation radius for the ACDMi (solid line) 
and ACDM2 (dashed line) cluster samples. If the perco- 
lation radius is significantly smaller than the mean ran- 
dom intercluster separation then the number of random 
supercluster is null and thus V ~ 1. As the percolation 
radius increases then the number of random superclus- 
ters increases as well, but with a slower rate than the 
corresponding number in a clustered distribution. How- 
ever, there is a trade-off which has to do with the fact 
that at extremely small percolation radii, although we 
ensure that V ~ 1, we get a very small number of su- 
perclusters. Therefore, the necessity of having adequate 
statistics leads to a larger percolation radius. Since in 
the present study we wish to investigate the morpholog- 
ical characteristics of superclusters and possible related 
environmental effects, we limit our analysis to super- 
clusters with more than 9 cluster members, for which 
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Figure 5. The M sup — a relation for ACDM rich superclus- 
ters (with 10 or more members). Left Panel. Based on rich 
clusters (ACDMi sample). Right Panel. Based on poorer 
clusters (ACDM2 sample). 



Figure 6. The frequency distribution of supcrclusters with 
4 or more members having their crossing times (along their 
minor axis) indicated along the x-axis. Left Panel. Based 
on rich clusters (ACDMi sample). Right Panel. Based on 
poorer clusters (ACDM2 sample). 



we can robustly determine their shapes (eg. Kolokotro- 
nis et al. 2002). We then find, for the percolation radii 
determined by eq.(l), ie., R pr ~ 18 and 11 Mpc 
for the ACDMi and the ACDM2 samples respectively, 
that the significance (which can be read from Fig. 4) is 
V > 0.95. 

An interesting property of our detected ACDM su- 
perclusters is a correlation between the supercluster 
mass, defined as the sum of the cluster members mass, 
and the cluster velocity dispersion within their parent 
superclusters, of the form M sup oc a a (sec Fig. 5). For 
superclusters with 4 or more members we find that the 
linear correlation coefficient in log-log space is ~ 0.4 
with the probability of random correlation being ex- 
tremely small (< 10~ 7 ) for both ACDMi and ACDM 2 
cluster samples. The corresponding slopes of the M — a 
relation is ~ 0.43 ±0.08 and ~ 0.56 ±0.06, respectively. 
Although the supercluster crossing time along their ma- 
jor axis is much larger than the age of the universe, t u , 
and thus superclusters have not yet turned around (see 
also Gramann & Suhhonenko 2002), such a correlation 
imply that they could be "locally" bound (eg. Barmby 
& Huchra 1998; Small et al. 1998). In fact, many super- 
clusters have crossing times along their minor axis less 
than t u - This is shown in Fig. 6 for the ACDMi (left 
panel) and ACDM2 (right panel) cluster samples. 

Furthermore, it is interesting that supercluster 
scale peculiar velocities can produce distinct features 
on the CMB and thus provide the means for detection 
of superclusters at cosmological distances (eg. Diaferio, 
Sunyaev & Nusser 2000). 



3.2 Application of the Shapefinder Statistic 

To identify the characteristic morphological features of 
the large scale structures, we will use a differential ge- 
ometry approach, called "shapefinders" and introduced 
by Sahni et al. (1998). 

The supercluster shapes are defined by fitting el- 
lipsoids to the distribution of cluster members. In 
Kolokotronis et al. (2002) we investigated the minimum 
number of points that are necessary in order to deter- 
mine unambiguously the true shape of a structure traced 



by a point process. Here we estimate shapes for those 
superclusters that consist of at least 10 cluster members. 

In the case of clusters we can define the shape di- 
rectly from the distribution of dark matter or gas parti- 
cles. In both cases we use the moments of inertia method 
(eg. Carter & Metcalfe 1980; Plionis, Barrow & Frenk 
1991). Summing up over the DM or gas particles (for the 
case of clusters) or the cluster members (for the case of 
supercluster) we define the inertia tensor: 



'a/3 = 



E 



(i) 



(2) 



which after diagonalisation provides the eigenvalues ai , 
02, and 03 and the corresponding eigenvectors which 
point into the direction of the principal axes. Since the 
eigenvalues can be considered as the three principal axes 
of the triaxial ellipsoid with the volume V = 4^ a 10203 
we can now relate the shape of the configuration to the 
shape of that ellipsoid. 

Below we present the main features of the shape 
statistic, following the notations of Sahni et al. (1998) 
(for a first application to astronomical data see Basi- 
lakos et al. 2001). To characterize the shape of any ob- 
ject these authors introduced a set of three quantities, 
defined by Hi = VS' 1 , H2 = SC' 1 and Hz = C, where 
V is the volume, S the surface and C the integrated 
mean curvature of the given object. These quantities 
have dimensions of length and are normalized to give 
Hi = R for a sphere of radius R. With those quan- 
tities one can define two dimensionless shapefinders: 
K=(K 1 ,K 2 ), where 



and 



H2 — Hi 
H2 ± Hi 



K 2 = 



Hz — H2 



(3) 



(4) 



Hz ± Ti.2 

The above technique characterizes the shapes of 
topologically non-trivial cosmic objects according to the 
following classification: 

(i) Pancakes if K 1 /K 2 > 1 

(ii) Filaments if Ki / K 2 < 1 
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(iii) Ribbons for {JKi,Ki) — (a, a) with a < 1 and 
thus structures K1/K2 — 1 and 

(iv) Spheres if cti = 012 — as and thus Ki = K2 = 

(v) ideal filament (one-dimensional objects) for Hi — 
H2 « Hz 

(vi) ideal pancake (two-dimensional objects) for 
Hi << H2 — H3 

(vii) ideal ribbon for Hi « H2 « Hi 

Ideal filaments (0,1), pancakes (1,0), ribbon structures 
(1,1) and spheres (0,0) are represented by the four ver- 
tices of the shape plane in the form of the shape vector 
K = (Ki, K2), whose amplitude and direction deter- 
mines the morphology of an arbitrary 3D surface. Fi- 
nally, for the quasi-spherical objects the -Ki,2 are small 
and therefore their ratio (K1/K2) measures the devi- 
ation from perfect spherical shapes. Note that below 
we will consider as spherical all simulated objects that 
have Ki, K2 < 0.005, the axis ratios of which lie in the 
range: 1 < ai/a 3 < 1.20 with (ai/oa) ~ 1.10 ± 0.06 
and 1 < ai/a 2 < 1.25 with (ai/a 3 ) ~ 1.15 ± 0.05. 

3.2.1 The Shape of DM and Gas Cluster Halos 

Using the above procedure we have estimated the shape 
of both the Dark Matter and the gas component of the 
ACDM clusters. We have performed this analysis for 
two minimal spanning tree linking parameters: for the 
nominal value of 0.17 and for 0.17/2" with n = 1 (prob- 
ing densities 8" times larger than in the former case and 
hence the inner dense cluster cores). In Fig. 7 we present 
the corresponding shape spectrum for the ACDMi (left 
panel) and the poorer ACDM2 cluster samples (right 
panel). The main panels show the results based on the 
Dark Matter cluster component while the inner panels 
to the corresponding gas component. The DM cluster 
component of both the ACDMi and ACDM2 cluster 
samples show a strong preference for filamentary-like 
(prolate-like) configurations with ~ 75% and 72% of 
the clusters having K1/K2 < 1 while only 1.7% and 
2.5% can be considered spherical. Note that in the main 
panels of Fig. 7 the shaded histogram corresponds to the 
larger linking length while the dashed line to the smaller 
one. The results based on the Dark Matter cluster com- 
ponent are very similar for both linking parameters in- 
dicating that both the inner cluster dark matter core 
and the cluster outskirts have a consistent shape. 

Regarding the gas cluster component the results 
are quite different. As can bee seen from the insert of 
Fig.7 the shape spectrum peaks at K1/K2 — 1 with 
little dispersion around this value. This result together 
with the quite small values of the Ki and K2 param- 
eters indicate that the gas distribution is much more 
spherical than the corresponding dark matter compo- 
nent, as expected. We find that the fraction of clusters 
having a spherical gas component is 53% and 64% for 
the ACDMi and ACDM2 cluster samples respectively, 
while for the smaller linking parameter (the dense clus- 
ter cores) the corresponding values are ~ 44% and 49%. 
This indicates that the gas distribution in the cluster 
cores is more flattened than in the cluster outskirts. Fur- 
thermore, the sphericity of the gas halos implies that the 




Figure 7. The cluster halo shape spectrum: Left Panel. 
Rich clusters (ACDMi sample), Right Panel. Poorer clus- 
ters (ACDM2 sample). The shaded histograms correspond 
to a linking parameter of 0.17 (p/pmcan — 330) while dashed 
lines to the denser (by a factor of ~ 8) core cluster region 
(linking parameter of 0.085). The insert shows the corre- 
sponding results for the gas cluster component. 




Figure 8. Gas halo sphericity versus the misalignment an- 
gle between the major axis of the DM and gas cluster halo 
components for the main linking parameter used (0.17). 
Left Panel. Rich clusters (ACDMi sample), Right Panel. 
Poorer clusters (ACDM2 sample). 



position angles orientations of gas halos are ill-defined 
and thus should be used in alignment studies with great 
caution. To visualize this effect we present in Fig. 8 the 
misalignment angle between the major axes of the DM 
and gas mass halos as a function of the ratio of the 
gas halo minor to major axis. Values of the latter near 
unity imply spherical gas halo shapes. It is evident that 
indeed the large misalignment angles are correlated with 
the gas halo sphericity. 

Finally, in order to test possible resolution effects 
on the halo shapes, we present in Fig. 9 a compar- 
ison of the DM and gas halo axis ratio distribution 
{M DM > 6 x 10 13 M Q ) in the 512 3 simulation with 
the same quantities obtained from another run of the 
same simulation with 8 times coarser mass resolution 
(2 x 256 3 gas and DM particles, see §2.1 for details) 
It is evident that the DM axis ratio distributions (left 
panel) are indistinguishable while there is a small effect 
for the case of the gas halos (right panel) , in the direc- 
tion of having artificially lower gas halo sphericity in the 
lower resolution simulation. 
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Figure 9. Comparison of the DM and Gas halo sphericity 
as a function of simulation resolution. In the left panel we 
also plot, as a continuous line, the Jing & Suto (2002) model 
which fits well the DM halo axis ratio distribution. 



0.8 



io.6 - 




15 16 17 18 19 20 10 12 14 

R pr (h M Mpc) R pr (tr> Mpc) 

Figure 11. The fraction of filamentary superclusters as a 
function of percolation radius used to define the superclus- 
ters. The arrow indicates the "optimal" percolation radius 
given by eq.(l). Left Panel. Rich clusters (ACDMi sam- 
ple), Right Panel. Poorer clusters (ACDM2 sample). 
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Figure 10. Comparison of the shape spectrum of superclus- 
ters containing 10 or more clusters (dashed-line histogram) 
and of the DM clusters for the two subsamples (left and right 
figure) . 



3.2.2 The Shape of Superclusters 

In Fig. 10 we show the shape spectrum of the detected 
superclusters (with 10 or more members - dashed- 
line histograms) and we compare with the correspond- 
ing spectrum of their cluster members (for a linking 
parameter of 0.17). It is obvious that the dominant 
shape of cosmic structures, being clusters or superclus- 
ters, is that of prolate structures (filaments), in agree- 
ment with previous theoretical and observational stud- 
ies (eg. Plionis et al. 1991; Basilakos, Plionis & Mad- 
dox 2000; Sathyaprakash et al. 1998; Basilakos et al. 
2001; Kolokotronis et al. 2002; Basilakos 2003; Einasto 
et al. 2003). Furthermore, it is clear that superclusters 
are more flattened structures than clusters themselves, 
populating a larger region in the shape-spectrum. 

We have verified the dominance of filamentary su- 
percluster shapes for the "optimal" percolation radius, 
given by eq. (1), although it appears that superclus- 
ters have a slightly lower filamentariness with respect to 
their cluster members. The relevant fractions are: fila- 
ments (~ 64%), pancakes (~ 32%) and ribbons (~ 4%). 
Note that there are no spherical superclusters. 

We have further investigated the dependence of the 
supercluster shape on the value of the percolation ra- 



dius and we present our results in Fig. 11 where we 
plot the fraction of superclusters having a filament-like 
shape (K1/K2 < 1). It is evident that for the richer 
cluster sample (left panel of Fig. 11), the filamentari- 
ness is a dominant feature only around the "optimal" 
percolation radius (indicated by the arrow). Similarly, 
for the poorer cluster sample (right panel) the "opti- 
mal" percolation radius coincides with a local maxima 
in the filamentary fraction. However, in this cluster sam- 
ple, the filamentariness appears to be a generic feature, 
almost independent of the percolation radius used. 



4 THE ENVIRONMENTAL TRENDS OF 
THE LARGE SCALE NETWORK 

4.1 Correlations of vector alignments 

As already mentioned in the introduction, an interest- 
ing question regarding environmental effects on large 
scales is whether clusters are aligned with the orienta- 
tion of their parent supercluster. To study the possible 
supercluster - clusters alignment we attach to each cos- 
mic structure (being either supercluster or cluster) the 
orientation as three unit vectors along the three main 
axes. 

The distance vector between the supercluster center 
of mass and the cluster member is r, and the normalized 
direction is f = r/r. Based on the notations of Beisbart 
et al. (2002) we consider the following alignment esti- 
mator (see also Stoyan & Stoyan 1994; Faltenbacher et 
al. 2002): 

A(r) = (|ei-e k |)(r) . (5) 

Therefore, for each supercluster I, the factor A(r) de- 
scribes the direct alignment of the corresponding vector 
ei with the vectors ek of the clusters within the par- 
ent supercluster. Note that A(r) is proportional to the 
cosine of the angle between ei and e k . The case of a 
random alignment signal implies A(r) = 0.5. 
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Table 2. Correlation analysis results for the different pairs x, y using the cluster DM component. The straight lines in Fig. 12 
represent the least-square fit to the data and have the general form y = Ax + ft. Also, 1Z is the correlation coefficient with the 
corresponding probability (Vr) of random correlation. Note that M sup has units of 10 14 ?t — 1 Mq. 



Sample 


Pair x, y 


A 


P 


n 


Vr 


ACDMi 


K X JK 2 , Aiax) 


-0.04±0.02 


0.63±0.02 


-0.39 


7.9xl0~ 2 




K1/K2, M aup 


9.3±4.2 


31.2±4.4 


0.45 


3.8xl0~ 2 


ACDM 2 


K1/K2, A(ai) 


-0.05±0.02 


0.60±0.02 


-0.29 


3.6X10" 2 




K1/K2, M sup 






-0.11 


0.40 



0.B - M„, ai.34xl0 14 h-'lyf; 




0.4 



• M„ 16.2x10'%-^; 





Figure 12. The shape-alignment (upper panels) and shape- 
mass (lower panels) correlations for the two supercluster sam- 
ples analyzed and at a = 0.8ai. The solid lines are the best 
least-square fits to the data (see Table 2). 



4.2 Shape- Alignment Correlation 

In Fig. 12 (top panels) we present the alignment A(ot\) 
estimator as a function of the supercluster shape, given 
by the ratio K1/K2. There is an obvious superclus- 
ter shape-alignment relation, seen in both cluster cata- 
logues (ACDMi and ACDM2), with supercluster align- 
ment increasing with supercluster filamentariness. Fur- 
thermore, only for the high cluster mass sample, there 
is also a significant correlation between the supercluster 
mass and shape, with mass decreasing with increasing 
filamentariness. Probably, this is to be expected since 
pancake-like structures, being two-dimensional should 
be more massive than the one-dimensional filaments. 
However this is not observed in the case of ACDM2. 

Table 2 summarizes the quantitative correlation re- 
sults for all tests using the cluster DM halo orientations. 
Probably due to the small number of superclusters (es- 
pecially in the higher cluster mass sample) the signif- 
icance of the correlations is not very high. In order to 
test whether the signal could be affected by a few cluster 
outliers at the supercluster edges, we have derived the 
correlation signal and its significance as a function of the 
parametrized distance of the member clusters from the 




Figure 13. The shape-alignment correlation coefficient (up- 
per panel) and random probability (lower panel) as a func- 
tion of supercluster-centric distance. 



supercluster center of mass (ie., as a function of r/ai, 
where r is the physical size of the supercluster and «i 
is its major axis). In Fig. 13 we present our results from 
which it is evident that the signal is relatively stable, 
although it appears to be higher and more significant 
around r/a\ ~ 0.8, indicating that indeed at larger dis- 
tances, some cluster outliers do reduce the correlation 
signal. 

We have also used the gas halo orientations to re- 
peat the same analysis but as expected from the results 
presented in section 3.2.1 (see Fig. 8) we do not find 
a significant signal for the ACDMi sample. However, 
for the poorer ACDM2 sample we have found a corre- 
lation signal (1Z = —0.25) with a random probability of 
V ~ 0.07. 

Furthermore, we have seen in Fig. 11 that the frac- 
tion of filamentary superclusters depend on the perco- 
lation radius used and thus if, as suggested by the cor- 
relations in Table 2 and Fig. 12, the alignment signal is 
higher in the filamentary superclusters we may expect a 
drop of the signal when the distribution of superclusters 
is dominated by pancake-like morphologies. Therefore, 
in order to investigate the above as well as the robust- 
ness of the correlation signal we re-assigned cluster DM 
halos to superclusters using a range of percolation radii, 
around the "optimal" value given by eq.(l). In Fig. 14 
we present both the value of the shape-alignment corre- 
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Figure 14. The shape-alignment correlation coefficient (up- 
per panel) and random probability (lower panel) as a func- 
tion of percolation radius used to define the superclusters. 
The arrow indicates the values corresponding to the "opti- 
mal" percolation radius. The graphs on the right and left cor- 
respond to the lower (ACDM2) and higher (ACDMi) mass 
samples. 



lation coefficient and its random probability as a func- 
tion of percolation radius used to define the superclus- 
ters. The resulting correlation coefficients is higher and 
more significant very near the "optimal" percolation ra- 
dius, defined by eq.(l), which also coincides with peaks 
in the distribution of supercluster filamentary fraction 
(Fig.ll). 

The increase of the supercluster filamentariness 
with alignment, is to be expected according to the hier- 
archical clustering scenario, where gas and galaxies flow 
into denser regions along anisotropic directions, defined 
by the large-scale structure (cf. West 1994; Ostriker & 
Cen 1996 and references therein). As a result of these 
inflows we expect clusters to be aligned with their neigh- 
bours, especially when both are members of the same 
filament. For low-f2 m cosmologies one expects that ma- 
jor merging and anisotropic accretion of matter along 
filaments will have stopped long ago. Thus gravitational 
violent relaxation would tend to isotropize the cluster 
phase-space configuration, more so in the recent times. 
However, it has been shown that clusters retain memory 
of their initial anisotropic configuration from were they 
accreted the main bulk of the matter that they contain 
(eg. van Haarlem & van de Weygaert 1993). 




Figure 15. The cluster velocity -supercluster alignment 
[t/(r)] and cluster major axis-supercluster alignment [A(r)] 
correlation for the ACDMi superclusters. 



where Vk is velocity vector with |vk| = 1, which corre- 
sponds to each host cluster while ei is the eigenvector 
of the major (r = cti) supercluster axis. 

We find that shape - velocity alignment correlations 
does not exist at any significant level. However, for the 
case of high mass clusters (ACDMi sample) there is a 
relatively significant correlation between A(r) and U(r), 
indicating that the alignment between supercluster and 
cluster member major axes is accompanied with a sim- 
ilar alignment between supercluster major axis and the 
direction of the cluster member peculiar velocity (see 
Fig. 15). 



5 CONCLUSIONS 

We have studied the morphological features of ACDM 
cluster dark matter and baryonic gas halos as well as 
superclusters using a large volume (V = 500 3 h~ 3 
Mpc 3 ) N-body+GADGET simulation with unprece- 
dented mass resolution. The measure of the struc- 
ture geometry has shown that prolate-like (filamentary) 
shapes dominates over pancakes, in agreement with 
other recent large-scale structure studies. We have also 
presented evidence that there is a specific link between 
the orientation of cluster dark matter halos and that of 
the large scale network. Cluster size halos appear to be 
aligned with their parent supercluster major axis, more 
so if the supercluster is filamentary-like. For the richest 
cluster halos we have also found a correlation between 
the halo peculiar velocity - supercluster alignment and 
the cluster major axis-supercluster alignment signals. 



4.3 Supercluster - cluster velocity correlations 

In this section we investigate possible correlations be- 
tween the supercluster shapes and the alignment be- 
tween supercluster major axis and their member cluster 
velocity field. Using the notations of section 4.1 we can 
define the following cross-correlation vector estimator: 

t/(r) = (|ei-v k |)(r) (6) 
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